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5_( Abstract 

The Regge- Wheeler equation describes Eixial perturbations of Schwarzschild metric in linear ap- 
proximation. Teukolsky Master Equation describes perturbations of Kerr metric in the same approxi- 
mation. We present here unified description of all classes of exact solutions to these equations in terms 
of the confluent Heun's functions. Special attention is paid to the polynomial solutions, which yield 
novel applications of Teukolsky Master Equation for description of relativistic jets and astrophysical 
explosions. 



1 Introduction 
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At present the study of different type of perturbations of the gravitational field of black holes (BH), 
^y-^ neutron stars (NS) and other compact astrophysical objects is a very active field for analytical, numerical, 

^ experimental and astrophysical research. Ongoing and nearest future experiments, based on perturbative 

f — and/or numerical analysis of relativistic gravitational dynamics, are expected to provide critical tests of 

the existing theories of gravity [T]. 

In the last five years the sensitivity of the operating detectors for gravitational waves LIGO, VIRGO, 
GEO, TAMMA has been improving at a formidable rate and one may expect the first direct observation 
of gravitational waves in the nearest future. Note that the existing observational results, collected in the 
last five years, give only limitations on the number of the BH-BH, BH-NS and NS-NS mergers. Since 
such mergers are still not observed, this number seems to be below the optimistic theoretical expectations, 
• • announced some eight years agcj^ The number of real mergers is believed to be consistent with the recent 

. theoretical and observational constraints [2 . The already started large projects like advanced LIGO and 

especially LISA, hopefully will bring into being the gravitational wave astronomy in the next decade. 
}-j Thus we are expecting to discover new fundamental physics. 

Another outstanding physical problem is presented by the gamma ray bursts (GRB) - the most 
powerful explosions in our universe after the Big Bang and the relativistic jets, related with them, as 
well as with other astrophysical objects. Due to recent developments of gamma ray astronomy in space 
missions SWIFT, Chandra, Huble Space Telescope, Spitzer, HETTE-2, BeppoSAX, and AGILE, together 
with ground observations by ESA and many other observatories, we already have very good observational 
data, which is still waiting for adequate theoretical explanation. The recently started Fermi/GLAST 
mission will give us more complete and precise data in the nearest future. Concerning the theoretical 
situation one has to stress that the existing theoretical models of GRB do not give a clear and acceptable 
explanation of the observational facts [3]. Moreover, there is a kind of crisis in this area, since the 
observational data seem to contradict the old models of central engine of long GRB. In addition, the 
presence of BH in short GRB was recently refuted by the existing detectors of gravitational waves jl] . 

New physical effects, due to the rotation of the gravitational field described by general relativity (GR), 
may play role in the supernova explosions. Mathematical tools for the study of the generation of waves 
with different spins during supernova explosions are needed, too. 

*Department of Theoretical Physics, University of Sofia, Boulevard 5 James Bourchier, Sofia 1164, Bulgaria, E- 
mail: fiziev@phys.uni-sofia.bg 

^Up to two detections of BH-BH mergers per year for LIGO I have been forecasted at that time. 
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A basic approach for theoretical study of the above phenomena is the perturbation theory of different 
background gravitational fields. At present we have two versions of this theory. One is based on direct 
perturbations of the space-time metric. The other one is based on perturbations of the scalars of Riemann 
curvature tensor, constructed using a specific tetrad in the space-time. In the present article we are 
considering the exact solutions of the perturbation equations in both versions. One of our goas is to 
compare the two sets of basic results. 

The well known Regge- Wheeler Equation (RWE) |2| 

+ {-dl + sUi) = (1.1) 

describes in linear approximation the axial perturbations of Einstein equations, using the Schwarzschild 
metric as a background. It plays an important role in modern perturbation treatment of Schwarzschild 
black hole (SBH) physics. Its study has a long history and significant achievements jHIH]. 

Usually one works only with bound ed so lutions of the corresponding angular equation [5j . Then the 



effective potential in the radial RWE (1.1 1 sUi{r) = (l — i) {^-^^i^ + has a specific dependence: 

Z (Z-|- 1 ) on an integer I > \s\. The area radius r > can be expressed explicitly as a function of the " tortoise" 
coordinate a; = r -I- Tg^,^ In (|r/rg^^ ^ 1|) using the Lambert-W function: r = LambertW(±e^~^) -I- 1. In the 
last formula the sign " stands for the SBH exterior r e (l,oo) and the sign " — " - for the SBH interior 
r £ (0, 1). Hereafter we are using units in which the Schwarzschild radius r^^i^ = 2M =1. In the potential 
sUi{r) the quantity s has values s = 0, ±1, ±2. The most important from astrophysical point of view are 
the cases s= ±1 and s= ±2, which describe electromagnetic and gravitational waves, correspondingly. 

The study of perturbations of the Schwarzschild metric was started in 1957 by Regge and Wheeler 
and was developed essentially by Chandrasekhar, Leaver and many others, see in [5l[ni[7|, in the review 
articles ^8^ and in the large amount of references therein. Especially, analytical study of the solutions was 
started in [7] and extended by different approximate methods, see in [S]. The exact analytical solutions 
of this problem have been found recently [9] . Using them one obtains the quasi-normal modes (QNM) of 
Schwarzschild metric under different boundary conditions in the most natural and straightforward way 
- solving numerically the corresponding boundary problems formulated in terms of the exact solutions. 
These solutions are useful both for a more deep understanding of the corresponding physical problems 
and for formulation of new ones. In the present article we give a more general consideration of the Regge- 
Wheeler perturbation theory, considering aZZ types of exact solutions of RWE. In particular, we pay special 
attention to the different types of polynomial solutions to RWE and recover new classes of such solutions. 
For this purpose we are to modify and extend to all cases the notions and notations, used in [9] for detailed 
description of the solutions. 

The ansatz s'^iit, f) = sRu>,i{i') with complex frequency lo ~ luu + icuj brings us to the stationary 
problem in the outer domain r > 1: 

dl + {io^ - sUi) sRu,i ^ 0. (1.2) 
Its exact solutions were described in |2] in terms of the confluent Hcun's functions HcunC(a, (3, 7, S, rj, z) 

m- 

In the exterior domain the variable r plays the role of 3D-space coordinate and the variable t - the 
role of the exterior time of a distant observer. Because of the change of signs of metric's eigenvalues 
At = gtt = 1 — 1/?' and = Qrr = ~1/(1 ~ l/^)i ™ the interior domain the former Schwarzschild time 
variable t plays the role of a radial space variable r^n = i e (—00,00) and the area radius r £ (1,0), 
plays the role of a time variable. As a result, the Regge- Wheeler "tortoise" coordinate x G (— oo,0) 
presents a specific time coordinate in the inner domain. Here for the study of the solutions of ( |1.1[ | 
it is useful to stretch the interior time interval to the standard one by further change of the variable: 
X ti„ = x—l/x € (—00,00). Thus we are placing the already existing singularities at natural places in 
the corresponding complex plane and facilitate significantly the numerical calculations ^ . 

These comments are important for the physical interpretation of the mathematical results. In partic- 



ular, the natural form of the interior solutions of (1.1) is: s*I'™;(''m, ^m) = e ^^^^"^ sRi^,i {r{ti„)), where 
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r{tin) = LambertW ( ~exp ( — \/ (%^) +1 ) ) + The dependence of this solution on the interior 



radial variable rj„ is simple. Its dependence on the interior time ij„ is governed by equation (1.2) with 
interior-time dependent potential sUi. Despite of this unusual feature of the solutions in the SBH interior, 
using this approach we obtain a basis of functions, which are suitable for the study of the corresponding 
linear perturbations. 

The negativity of the imaginary part lui = < of the frequency lo ensures linear stability of the 
solutions with respect to the future time direction t — > -|-c» in the exterior domain [Sllllj. This condition 
is not enough to guaranty the stability of the interior domain (See detailed analysis in terms of the exact 
solutions to RWE in 0.). 

Usually one is not considering the domain r < 0, despite of the fact that it is not excluded by the 
geometrical meaning of the area radius r. In this domain the variables t and r restore their original 
meaning. We have to point out that for the study of the analytical properties of the solutions of RWE 
one must consider complex values of the area radius r € ^ . 

The study of perturbations of rotating relativistic objects in Einstein GR was pioneered by Teukol- 
sky [12] making use of the famous Teukolsky Master Equation (TME). It describes the perturbations 
s'tlj{t,r,9,(p) of all physically interesting spin-weights s = 0, ±1/2, ±1, ±3/2, ±2 to the Kerr background 
metric in terms of Newman-Penrose scalars. The pairs of spin-weights s with opposite signs a = sign(s) = 
±1 correspond to two different perturbations with opposite helicity and spin \s\ — 0,1/2,1,3/2, or 2. 
Under proper boundary conditions for TME one obtains QNM for Kerr black holes (KBH). The various 
significant results and references may be found in [BJ |5] . 

We have to stress one general feature of the description of rotating relativistic objects. The Kerr 
metric describes exactly the vacuum solution of Einstein equations - KBH. The gravitational field outside 
the rotating compact matter objects differs form the one, described by Kerr solution in higher order 
multipole momenta, due to the corresponding specific matter distribution (mass distribution and current 
distribution) [13J . In the outer domain the contribution of the l-th multipole moment to the gravitational 
field of given object is of order ^(^)', I = 0,1, 2,.... Therefore in the outer domain one can fit the 
gravitational field of any rotating compact object by Kerr metric with proper parameters M and a. This 
fitting is exact for the first two multipole terms. The difference between the real metric and the Kerr one 
will appear only in the next terms, which decrease with the distance as r^'^^, I > 2. Hence, the effects 
of the real matter distribution inside the compact object are negligible at distances greater than several 
event-horizon-radii and one can use the Kerr metric outside the event horizon as a very good approximation 
to the gravitational field of rotating compact objects of different nature. This is the physical basis for our 
applications of solutions to the Teukolsky Master Equation for study of different real relativistic objects 
like massive compact objects, jets, supernova explosions, etc. 

On the other hand, the above argument shows that it is impossible to make a decision about the true 
nature of the observed massive dark compact objects like ABH, studying their non-perturbed gravitational 
field outside the event horizon and extracting only the parameters M and a from the observational data 
for the first two multipole momenta. The same values of the parameters M and a can be prescribed to 
rotating relativistic objects of different kinds. 

The contemporary observational data is not enough to decide among the possible alternative interpre- 
tations of massive dark compact objects in astrophysics (See for example the very recent review by Matt 
Visser in We need reliable methods, which are able to make difference between models, like BH, 

naked singularities, superspinars, gravastars, boson stars, soliton stars, quark stars, fuzz-balls, dark stars, 
etc, with the same mass . These methods have to give indisputable evidences for the real nature of the 
observed compact dark objects, currently named astrophysical black holes (ABH). 

Such methods, based on the study of spectra of perturbations of the gravitational field of the real 
astrophysical objects were proposed in [HIH]. Using them one can "see" directly the event horizon, if it 
exists, due to the unique boundary conditions on it. In |14j only approximate methods for study of these 
problems were discussed. The exact solutions to the RWE and TME, studied in the present article can 
help to extend the approach of [S] to the problems with rotating relativistic objects. 

The key feature of the TME is that in Boyer-Lindquist coordinates one can separate the variables using 



3 



the ansatz tjj{t,r,d,Lp) — e~*"*e™'^S'(6')i?(r), i.e. looking for solutions in a specific factorized form. From 
mathematical point of view the function %{t,r,6,(p) ~ e~*'^*e™'''S'(6')i?(r) actually defines a factorized 
kernel of integral representation for the solutions to TME: 



m— — oo 



basically introduced in the problem at hand for the first time in [T^] . The study of the QNM [8j showed 
that this kernel can be singular with respect to the variables E and r: it is proportional to Dirac d{E—Esp), 
where Egp belongs to some specific for the given problem spectrum; with respect to variable r the kernel 
may be singular at infinity and at the horizons. In the existing literature only regular with respect to the 
variable 6 kernels are in use. In the present article we are starting the consideration both of regular and 
singular with respect to 9 kernels in the natural integral representation ( |1.3[ ) of solutions to the TME. 
The different type of kernels are to be used for solution of different boundary problems. Note that from 



physical point of view is important the regularity of the very solution s'>p{t,r,9,(j)) in equation (1.3 1. The 
kernels like 3C(t, r, 9, f) are an auxiliary mathematical objects. One is often forced to use singular kernels 
for the natural integral representations of the solution of physical problems. The regularity of the function 
sV'(ij 1^, (/>) with respect to the variable 9 depends on the choice of the amplitudes s^cj,_E,m- 

As a result of factorization of the kernel 3C(t, r,9,ip) a pair of two connected ordinary differential 
equations for its nontrivial factors sSuj,E.m{9) and sRu],E,m{f) arises - the Teukolsky angular equation 

(TAE) [niiii] 

'sin9^.S^,E^m{9)\+ sW^,E,m{9)sS^.E.Ad)=Q, (1.4a) 



sin ed9 \ d9 

sW^^E.rn{9) = + a^w^ cos^ 6* - 2saLucos9- (m^ + + 2tos cos 6*) / sin^ 9; (l-4b) 
and the Teukolsky radial equation (TRE) [12] 



A-^|- (^A^+i^ ,i?^,B,™(r)^ + ,K.,£,™(r) .i?c.,B,™(r) = 0, (1.5a) 
A A ar 

Here the azimuthal number m has arbitrary integer values m = 0, ±1, ±2, . . . , and A = — 2Mr + a^, 
K = aj(r^ + a^) — ma, L — E — s{s + 1) + a?uP — 2mauj — 4 isuj r. The standard real parameter a is related 
with the angular momentum of the Kerr metric, M is the Keplerian mass of the Kerr solution. The two 
complex parameters lu and E - the constants of the separation, are to be determined using the boundary 
conditions of the problem. 

As in the case of Schwarzschild background, the negativity of the imaginary part luj — "^(lu) < of 
the frequency lo — luj^ + iui may ensure linear stability of the solutions in the exterior domain of the Kerr 
metric with respect to the future time direction t +oo |12l I16j . In the interior domain the solutions to 
TME are not stable [17 . 

Much like in case of the Schwarzshild solution, in the interior of the Kerr metric - between the zeros 
r±: Q<r_<r^oi the function A, two of the eigenvalues: and of the metric in Boyer-Lindquist 
coordinates simultaneously change their signs. Indeed, one pair of eigenvalues are \g — ggg — r^ + a^ cos^ 9 
and \r = grr = {r'^ + a^cos^6')/A. The second pair of eigenvalues are the roots At, A^ of the equation 
•^^ " (.9tt + 9(i>4>)^ + 9tt94>4> — Qt^- Their product equals AjA^ = — Asin^ 9. The last expression, together 
with the form of proves the simultaneous change of the signs of the two eigenvalues Af,Ar, when 
the variable r crosses the horizons r±, since the determinant of the metric g ^ — (r^ + cos'^ 6*)^ sin^ 6* 
does not vanish there. As a result, between the two horizons r± the variable tin = x € (—00,00) plays 
the role of the interior time and the variable ri„ = t is the interior radial variable. We are using the 
following Kerr-metric-tortoise-coordinate: x = r + a+ In |r/r_|_ — 1| — a_ In \r/r^ — 1| G (—00, 00), where 
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a± = r^~^r It is a straightforward generalization of the proposed in [T^] tortoise variable for the 
exterior domain r G (r+jOo). Since our expression is valid in the interior domains, too, the inverse 
function defines r = r(tin) when r G (r_, r_|_). In the second interior domain r < r- the variables r and t 
restore their original meaning. For a detailed analysis of the light cones in Kerr geometry see |18j . 

Despite of the essential progress both in the numerical study [19] of the solutions to the equations 

at present there is a number 
that RWE 



(1.4a I and (1.5a I and in the investigation of their analytical properties 



of basic questions remaining unanswered. For example, it is well known for a long time 



1.2h, as well as TAE (1.4a I and TRE (1.5a I can be reduced to the confluent Heun differential equation 



H" 



+ 1 , 7 + 1 

z z — 1 



H' 



1 



H = Q. 



(1.6) 



Here the constants and v are related with the accepted in the notation HeunC(a, /3, 7, (5, t;, z) ones: 
a, (3, 7, 6, T] according to the equations 



S — ^ + V - 



(3 + J + 2 



a{(3+l) /3 + 7 + /37 
^= n M • 



(1.7a) 
(1.7b) 



To the best of our knowledge the exact analytical solutions of the angular and radial equations ( 1 .4a ) and 



(1.5al are still not described in the literature in terms of confluent Heun function HeunC(ck, /3, 7, (5, 77, z) 
- the unique particular solution of the (1.6 1, which is regular in the vicinity of the regular singular point 
z = and normalized by the requirement HeunC(a, /3, 7, S, 77, 0) = 1 (See for details [lOj.^ 

In the late 2006 a program for flUing this gap has been started as a natural extension of the articles 
[5]. The first results were quite stimulating j22j, but serious difficulties came across both in analytical 
and in numerical studies. This is because the theory of Heun's functions, as well as numerical tools for 
computer calculations with them are still not developed enough. The main purpose of the present series 
of articles is to report some of the basic results, obtained for detailed description of the exact solutions 
and to describe different boundary problems for RWE, TRE, TAE, and TME in terms of the conffuent 
Heun function, developing both the theory and computational tools for this function. 

We outline several applications of these solutions to well known physical problems like QNM of static 
and rotating BH, as well as applications to other astrophysical phenomena. In particular, we found novel 
simple mathematical description of relativistic jets, applying specific solutions of TME, which were not 
used up to now |23j . We have to stress that the relativistic jets are quite common phenomenon at very 
different physical scales |24]. Around very different astrophysical objects: small mass brown dwarfs, 
neutron stars, protostars, GRB, quazars, blazars, galactic centers, and even around galactic clusters have 
been observed jets of various size, corelated with the scales of these objects. Therefore, if possible, a 
common description and models of relativistic jets are highly desirable. 

According to [TU] the confluent Heun function HeunC(a, (3, 7, S, ?], z) reduces to a polynomial of degree 
iV > of the variable z, if and only if the following two conditions are satisfied: 



S /3 + 7 
a 2 



+ A^+1 = 0, 

AAr+i(/i) =0. 



(1.8a) 
(1.8b) 



a " 5-condition" , and the second one, i.e., the ( 1.8b[ ) 



Further on we call the first condition, i.e., the ( l.i 
- a "Ajv+i-condition". 

^ In the present article we are using the Maple-computer-package notations for Heun functions. Basically, these notations 
are borrowed from the two mile-stone articles on modern theory of Heun functions by Decarreau et al. in |10l . and at present 
seem to be most popular, since the Maple package is the only one for analytical and numerical work with Heun functions. 
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We represent the three-diagonal determinant Ajv+i(/^) in the following specific explicit form: 



^l-ql + 

Na fi-q2 + la 2(2+/3) 
(iV-l)a fi-q3+2a 





















(1.9) 





















^i-qN-i + {N-2)a {N-l){N-l+(3) 
2a n-qN + iN-l)a 

la 





N{N+(3) 
fJ'-QN+i+Na 



which turns to be useful for calculations. Here qn = (n — l)(n + + 7). 

On the other hand, the so called algebraically special solutions of RWE and TRE were discovered 
long time ago |25j . These are of a generalized polynomial type. According to the existing literature these 
solutions describe pure incoming or pure outgoing waves. To the best of our knowledge the algebraically 
special solutions are still not discussed in terms of Heun polynomials. 

As far as we know, applications of this class of solutions to real physical problems can not be found 
in the existing literature on gravitational physics. If such waves are spreading only in one direction, they 
seem to be most suitable for description of relativistic jets and supernova explosions. 

We have to stress one more argument in this direction. Very recently the algebraically special solutions 
of RWE and TME were proved to be relevant for the study of instabilities of different kind of some more 
or less "exotic" solutions to the Einstein equations Obviously, the physical manifestation of the 

instabilities of the mathematical solutions are the explosions of the corresponding objects. Therefore it 
seems natural to look for a description of astrophysical explosions in terms of polynomial solutions of 
TME, which are stable in the future and instable in the past. 

Our study of the polynomial solutions of TAE, which were completely ignored in the gravitational 
physics up to now, shows that these indeed may describe in a most natural way the coUimation of the 
jets and demonstrate reach variety of jet forms, very similar to the observed ones [HJUS]. We are going 
to analyze these problems on a correct mathematical basis using the confluent Heun functions. 

It is well known that the Kerr metric approaches the non-rotating Schwarzschild one in the limit a — > 0. 
One of our aims is to stress the analogy and the differences between the solutions of RWE and TME. These 
two equations describe the perturbations of Schwarzschild and Kerr metrics in terms of different quantities: 
RWE - in terms of direct perturbations of metric and TME - in terms of perturbation of Weyl scalars. 
Therefore in the limit a — s- there exist a smooth transition from perturbations of KBH to perturbations 
SBH in terms of Weyl scalars, but a simple transition from the solutions of TME to the solutions of RWE 
is not possible [6l. Nevertheless, the mathematical analogy between the corresponding solutions becomes 
quite transparent when the solutions are represented in terms of confluent Heun functions. 

Another point of ours is to collect at one place the known facts, together with the new results in the 
field and to describe all of them uniformly and in common terms and notations. This way we hope to 
obtain a more clear picture of the quite complicated present-day state of the arts in the perturbation 
theory under consideration and its possible further developments. 

The large amount of the accomplished work and the number of the obtained results forces us to publish 
them in a series of articles, starting with this one. 
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Exact Solutions to the Regge- Wheeler Equation in Terms of 
the Confluent Heun's functions 



In area-radius- variable r the RWE (1.2 1 reads: 



1 

1 - r 



dr 



(2.1) 



+ (/a + i)-.^ + i) [1 + ^) 



- 1 




0. 



Two of them: 
oo is an irregular 



The anzatz sRi^,i{'f) = r*+^(r — l)*"e'"''iJ(r) reduces it to the confluent Heun equation ( |1.6[ ) with the 
following specific parameters: a = 2iuj,P — 2s,j = 2iw,6 = 2uj^,ri = — l{l + 1) [5]. This anzatz was 
used for the first time for analytical and numerical study of the solutions of RWE by Leaver in [7j without 
reference to Heun equation and Heun functions. 

The equation (1.2 1 has three singular points in the whole complex plane p, [TH]. 
r = and r = 1 are regular and can be treated on equal footing. The third one: r 
singular point. 

Note that, after all, the horizon r = 1(= 2M) turns out to be a singular point for the perturbations of 
SBH, despite of the fact that it is not a singularity of the algebraic invariants of the Riemann's curvature 
tensor 'Jlijki- The algebraic invariants are able to indicate only the curvature singularity at the point 
r = 0. The horizon does not define a singular surface of these invariants in the Schwarzschild space-time 
manifold M*^^'"^) and usually is treated as a pure coordinate singularity of the metric in Schwarzshild 
coordinates. Since the algebraic invariants do not fix the whole geometry of the space-time manifold, their 
consideration is not enough to recover all space-time properties. For this purpose one must consider a 
large enough number of high-order-differential-invariants of the Riemann's tensor Hiju [27| . For example, 
the differential invariant of first order 



DI 



(V In r) 



1 

36 



(Vln 



1 



(2.2) 



can be used as a coordinate-independent indicator both of the event horizon r = 2M{— 1) of SBH and 
of the singularity at r = 0. Obviously, both of the last two values of the area variable r define non- 
coordinate geometrical objects in the Schwarzschild space-time manifold M^^''^' and at the same time - 
alike singularities of equation (2.11. It becomes clear that for the study of small deviations from the 



background metric it is necessary to consider differential invariants of the Riemann's tensor, relevant to 
corresponding variations of geometry. 

Using the confluent Heun function one can write down 16 local Frobenius type solutions of the RWE 



(2.1) 



7->± 



(r)/r : 



HeunC(craQ;^ , ap(i^ , cr^7^ , (5^ , 77^ , ). 



(2.3) 



Here 



= 2iuo^ 13 . = 2s, 7, — 2iuj, S. ~ 2a>^, 77 



a _ = -2iw, /3_ = 2ia; , 7_ = 2s, (5_ = -2lj2 , 77_ =2lj2 + - / (/ -I- 1 ) 



±1,CT^ 



±1 



According to equations (2.3) and (2.4d) the behavior of the solutions sR^i, 



(2.4a) 
(2.4b) 
(2.4c) 
(2.4d) 

(r) around the 



corresponding singular points z± = is deflned by the dominant factor z^^^"^*/^. All other factors in (2.31 
are regular around these points. In contrast, the solutions sR^ i „ „ ^ (r) are in general singular around 
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the corresponding singular points z± — 1, i.e., around r = - for i 



(r) and around r ~ 



for ,r: 



(r). This explains the meaning of the upper ± signs in the notation 



Only two of the sixteen solutions (2.3 1 are linearly independent. For different purposes one can use 



different pairs of independent local solutions. 

In addition, using the asymptotic expansion of the confluent Heun function [lOJ we obtain two asymp- 
totic solutions of Tome type, i.e. the local solutions around the irregular singular point \r\ = oo of the 
RWE in the complex plane 



jDltoo 



( cj(r-l-lnr) ^ ^ ^ 



±iy 



flo 



(2.5) 



For the coefficients aj — aj,u.i,iya,<yi3,a^ one has a recurrence relation |10j which shows that they increase 
together with the integer j. Hence, the asymptotic series (2.5 1 is a divergent one. 



3 Polynomial Solutions to the RWE 

We have to remaind the reader that according to the accepted in [10] terminology " polynomial" are called 
not only the cases in which the confluent Heun function is indeed a polynomial. The same terminology 
is currently in use for all solutions which have a form of finite Taylor series expansion multiplied by 
elementary functions. We shall apply this terminology to the solutions ( |2.3| of the RWE and TME, too. 
This way we obtain a more reach class of quasi-polynomial solutions |10l I25j . 



The 5-condition (1.8a I gives 





r s 1 






cr/3 < 









+ N + 1= 0, N>0; 
and yields the following two classes of polynomial solutions to the RWE. 



(3.1) 



3.0.1 First Class of Polynomial Solutions to the RWE: 



For the solutions , 
CT/3 = — (T where a 



R^, 1 „ „ „ {f) with a„ — 
— sign(s). The same result 



cr^, and Op 



±1 equation (3.11 gives 



+ 1 = \s\ we obtain for the solutions 



.r: 



sR: 



.La, 



(r) 



R+ 



, if 
(r) 
(r) 



when (Tq = cr^ and = —a. For brevity we will use the short notations 

and sR~, 1 „ (r) — sR . i „ „ „ (r)- 

Applying the described in the Introduction general construction, we see that if in addition the A^v+i 
condition 



(|1.8b|) is fulfilled in the form Aj^|(^) 



0, then the confluent Heun function in (2.3) reduces 



to a polynomial of degree |s| — 1 > 0. The A^v+i-condition gives a finite number \s\ different eigenvalues 
^J'■k=l,....\s\ of the parameter fi for the problem at hand. 

1. In the case of RWE in general \s\ = 0, 1, 2, but the requirement sN = |s| — 1 > is satisfied only 
for \s\ = 1,2. Hence, first class polynomial solutions of RWE may exist only for electromagnetic (|s| = 1) 
and gravitational (|s| = 2) waves and there are no such solutions of scalar (|s| = 0) nature. 

2. The explicit form of the Ajv+i-condition, derived from (1.8b I for electromagnetic waves, is A^(/i) — 
jjL — Q. Then one obtains l{l -I- 1) = 0, making use of the general relation between /i and r] - (1.7b I and the 



expressions for r]± in (2.4a I and in (2.4b I. Hence, polynomial solutions of the first class for radial RWE 
with \s\ = 1 and integer I > 1 actually do not exist. The regularity requirement, usually posed on the 
solutions of the angular RWE |3-|H], is obviously too restrictive for existence of first class electromagnetic 
waves with polynomial dependence on the radial variable r. 



^The notation ±oo in (|2.5 
singular point |r| = oo, i.e 



notation z±, since z± 



indicates the two different directions on the real r-axes in Cr 
oo denotes the limit r +oo and — oo denotes the limit r 
itoo when r — ► oo. 



for approaching the irregular 
-> — oo. It is consistent with 
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3. The explicit form of the AAr+i-condition, derived from (1.8b I for gravitational waves, is 



^2 (W = M -2(1+ ^ — — ) /i - cr„a± (1 + cr/3/3±) = 0. 



(3.2) 



Using the roots of this equations, the (1.7b), the expressions for ri± in (2.4al and in (2.4b), one obtains 
the following spectrum of the polynomial solutions sR^ i „ (jY 



2^1,, 



cT^-{l-l)l{l + l){l + 2), ; = 2,3,4,.. 
6 



(3.3) 



This spectrum describes stable in the future [t —>■ +oo) gravitational waves when CTq = — 1, and stable in 
the past (t — oo) gravitational waves, if cTq = +1. 

Obviously, using the properties of the confluent Heun function this way we have re-derived the well 
known algebraically special solutions to RWE [25 , which describe a specific class of one-way gravitational 



3.0.2 Second Class Polynomial Solutions to RWE: 



Let us impose the (5-condition on the other solutions in (2.3). Then because of the asymmetry between 



the parameters f3 and 7 in (1.7b I we have two different subcases 



A. In the first subcase we consider the solutions ci? , , „ „ „ (r) with (t„ = —a-,. We denote these 



solutions as sR^ i a ~ sR^ i -a af, <T ('')■ them the (1.8al gives an infinite series of pure imaginary 
equidistant spectrum of frequencies 



-{N+l + aps), iV-0,1,2, 



(3.4) 



Making use of (1.7b I we obtain for the parameters in (2.4a) 

sl^NEa,a,=E+{l + <ypsf + {l+2<7ps)N, TV = 0, 1, 2, . . . , 

with E = l{l + 1), where / = \s\, \s\ -I- 1, |s| -I- 2, . . . is an integer. 



(3.5) 



B. In the second subcase we consider the solutions ci? 



(r) with (Tq = —cTfj- We denote them as 

' I a — s-^'-Lj I -a/i cr^ a ' ^'^^ them the ( I.Sal gives an infinite series of pure imaginary equidistant 
spectrum of frequencies 



.r: 



= .r: 



s^Na,a,^^Pl.iN+^+^l^)^ iV = 0,l,2,, 



Making use of (1.7b I we obtain for the parameters in (2.4b I 

s^^NEa,a-,^ E + a^sN +{N +l + G^sf, iV = 0,1,2,. 



(3.6) 



(3.7) 



with E = 1(1 + 1), where I — |s|, \s\ -I- 1, |s| + 2, . . . is an integer. 

li E = l{l -1-1), we will have regular solutions to the angular equation - the spin-weighted spherical 
functions. It is easy to check that in this case the AAr+i-condition can not be satisfied for N = 2 and for 
> 4. In contrast, for = 1, 3 we have several new polynomial solutions to the TRE: 
1. The A /V4_i -condition is fulfilled for the solutions ci?^, (r) when = 1 and: 



a) s = 2, ? = -2, 1, aa = Tl, ct/? = -1, cr-. 



±1; 



b) s 



-2. I 



-2, 1, (Ja = Tl, (^(3 = 1, 0^7 = ±1- 



2. The AAr_|_i-condition is fulfilled for the solutions sR^ 

a) s=\,l^ —1,0, cFa = 1, cT/g = ±1, cr-v = —1; 

b) s = —1, / = —1,0, (Jq, — —1, (7/3 

c) s = 2, I = —2, 1, Ga — ^, <yp — ±1, cr 



V) when A'^ = 1 and: 



7 — 
±1, (T^, 
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d) s = —2, I = —2, 1, (Tq, — —1, ap — ±1, a-y — 1. 

Since these solutions do not satisfy the requirement ^ > |s| for the solutions to angular equation, they 
seem to be incidental. We include them for completeness of our list of polynomial solutions to the radial 
RWE. 

3. The A AT+i -condition is fulfilled for the solutions (r) when N = 3 and: 

' ^ ,t -CTq ,fT^ ,cr-y V / 

a.) s — 1, I = —2, 1, CTq = ±1, (Jp = —1, (T-y = =f1; 
b) s = —1, I = —2, 1, (Jq, = =f1, (t^ = 1, ct^ = ±1. 

The solutions with = 3 and I = —2 do not satisfy the additional condition I > \s\, but for / = 1 
we obtain two new algebraically special solutions to the RWE of electromagnetic type with = ±i3/2, 
which are not described in the known to us literature. 

The situation with the polynomial solutions to the RWE will be drastically changed if one rejects the 
regularity requirement on the solutions to the angular equation. As a result the condition E — 1(1 + 1) 
will be no more valid. Then in the case of polynomial solutions of second class one can use the (13. 4| and 



(3.5), or the (3.6 1 and (3.7 1 in the A7v+i-condition. Thus one obtains an algebraic equation of degree 



(N + 1) for the constant E with solutions E = ^E^ ^ ^ ^ ^ where n — 0, . . . , N 



The infinite series of polynomial solutions of second class n o- o-^ o- ('')' entering the expression 
s^N,n,a^,a0,a^{t,r) = e'*'^"'"'"" '"-3 '"^ ^^^^ (r) , dcscribcs non-oscillating exponentially decaying 

in time linear perturbations of SBH, when 5 ^s'^^o-^o- ^ < 0. Obviously, the infinite set of these solutions 
for A'^ = 0, 1, 2, . . . presents a complete basis for Laplace series expansion in time t for linear perturbations 
of more general form. The angular factor gS*^ ^ ^ ^ (9) of the corresponding kernel 

will be singular at least at one of the poles 9 — 0, 9 — tt, or simultaneously at both of them. This causes 
a strong anisotropy of the emission of the corresponding waves, i.e., a collimation of the emission around 
the poles. This kernel can be used for representation of this type of solution to the full RWE in the form 

{t,r,9,^). (3.9) 



The regularity of these solutions and the convergence of discrete sum (3.9 1 depends on the choice of the 

amplitudes AjM^n,7n,cr^,cri,,a^- 

One can find analogous finite set of polynomial solutions of the first class to the radial RWE, which 
correspond to singular angular part of the total linear perturbation, too. 

We will not discuss here these r-polynomial-0-singular solutions in detail, since it is more convenient 
to obtain the corresponding physical information in terms of perturbations of Weyl scalars using TME in 
the limit a — > 0. This will be done in the next Sections. 

4 Exact Solutions to the Teukolsky Radial Equation in Terms 
of the Confluent Heun's functions 



The explicit form of TRE 



CpRui E.m , f-, , \ f ^ , 1 A dRuj.E,m 

-f (1 + s) ' ' ' 



dr^ \r — r+ r — r- ) dr 

2 



(uj (a^ + r^) - arri\ / ^ 1 \ / \ 

\ u V- - + r + - 

(r — r^. j (r — r_ j \r — r — J \ / 



-E + s{s+l)-a^u>'^ + 2mauj + Aisujr\ -, ''^^ . (4.1) 

(r — r+)(r — r_) 
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shows that it has three singular points: r — r± and r ^ oo. In the present article we consider only the 
non-extremal Kerr metric with real r+ > r^ > 0. Then precisely as in the case of RWE, the first two 
are regular singular points (the Cauchy horizon r = r_ and the event horizon r = r_|_), and the third one 
(the physical infinity r = cx)) is irregular singular point. The symmetry of the (4.1 ) under the interchange 
r+ ^ r_ is obvious. 

As well as in the case of SBH, the algebraic invariants of Riemann tensor are not able to indicate the 
horizons of KBH and one usually considers them as a pure coordinate singularities of the metric in Boyer- 
Lindquist coordinates. In contrast, the circle r = 0,9 = tt/2 is a singularity of the algebraic invariants of 
Riemann tensor [6]. 

It is not too hard to find differential invariants of the Riemann tensor, which are able to distin- 
guish both horizons r — r± and the ergosphere gu = 0. Indeed, let us consider the following alge- 
braic invariants of Weyl tensor WijkV- Ii — j^'V^ijki^^''^^ ~ the density of the Euler characteristic class, 
and I2 — jgy^ijki*^^ ^ the density of the Chern-Pontryagin characteristic class [28]. Let us put 

/ \l/6 / \l/6 

[h - ihY'^ = A = |A|exp(i?A). Then r = i^A cos(?/'/6) and p = [J^A cos(?/'/6)-i are obviously 

invariants of Weyl tensor - non-algebraic and non-differential ones. In Boyer-Lindquist coordinates one 

2 

obtains p — r + ^ cos 9 and gu = 1 — 2M/p. The differential invariants of first order 



DIi = -(Vlnr)^ = — 
rp 



1 - 



2M 



DI2 = (Vlnp)^ - (Vlnr) 



1 1 



a 

-2 

2M 
P 



(4.2a) 
(4.2b) 



indicate the two KBH horizons, the ergo-sphere and some other geometrical objects in Kerr space-time. 



The very horizons are singularities of the same kind in the equation (4.1 1. The relations in (4.2) generalize 
the equation (2.2 1, which is a limiting case of equation (4.2a I for a — > 0, when p r. In this limit the 



invariant in equation (4.2b) becomes trivial: DI2 0. Hence, at this point we have a complete analogy 



between RWE and TRE. 

The analytical study of the solutions to the TRE and TAE was started in [7 and continued by 
different approximate methods (see in [HJIin]) without utilizing of Hcun functions. Using the confluent 
Heun function one can write down 16 exact local Frobenius type solutions to the TRE (4.1 1 in the form: 



{r;r+,r^)A'/^^e° 



± T 



which is very similar to the form of the equation (2.3). NowJ^ 



HeunC(CTaa^ , , 0-^7^ , (5^ , 77^ , ), (4.3) 



a^^sau,,E,ni{r+,r^)=2iuj{r+ - r_)= ipuj/fla, 

/3+ =s/3c,,E,m{'^+,r-)= s + 2i{m-uj/n_) /p, 

7+ =sluj.E,m{r+,r^)^ s - 2i{m - lo/VL+) /p, 

<5+ =s5uj,E,m{r+,r^)^ a+ (s - iuj{r+ -|- r_))= a+ (s - iuj/Vtg) , 



'n+=s'nu,,E,m{r+,r^)-- 
z+^z+{r]r+,r^)^- 



-E 



2m^fl^ 



(2mfla 



p^m 



r^ — r 

r-\- — r^ 



z_ = z_(r; r_|-, ?'_) = — ; z_|_ + z_ = l, 



UJ f2 + 
'i 

-A 
(r+ — r_)^ 



(4.4a) 
(4.4b) 
(4.4c) 
(4.4d) 

(4.4e) 
(4.4f) 



The discrete parameters aa,crf3,(Jj are the same as in (2.3). In equations (4.4a)-(4.4e I we are using 
the following quantities: the angular velocity of the event horizon f2-|_ — ^/rlT/r^ / (r+ + r_), the angular 



^Note that the notations z± in the ( |4.4f| are consistent with the corresponding ones for the case of RWE - ( |2.4d[ l. They 
are based on the hmits z± — > ±00 for r — > 00. Their relation with the notations of the parameters of KBH r± is illustrated 
by the equations 2±(r:p; r+, r_) = 0. The labels it in the notation — \4.3\ are related to the labels of their arguments 
z±, not with the labels of the parameters r±. 
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velocity of the Cauchi horizon ri_ = yJr^^^JrZ j [r + + r -) , the arithmetically-averaged angular velocity 
= (ri_|_ + ri_) /2 = l/(2a), the geometrically-averaged angular velocity fig = ^JT^JlZ ~ 1/(2M), 
and the new dimensionless parameter p = yjr+jr^ — \fr^^Jr+ — \/^-/^+ ~ G (0, oo). Note 

that the inverse relation r± = \/flzf/il±/ + ri-) permits us to replace r± with il± wherever it is 
necessary, thus making transparent the duality of the parameters r± and fl±, as well as the behavior of 



-P, ^a., 



the above quantities under interchange of the two horizons: 
- invariant. 

The parameters a_, P_,^_,S_,ri_ can be obtained by interchanging the places of the two horizons: 
r+ ^ r_ in (4.4al - (4.4el. This procedure may be substantiated using the known properties of the 
confluent Heun function under changes of parameters |10| . One can check directly that in this way we 
obtain indeed solutions of (4.1 1. 

.a„.aAr;r+,r_) 



According to equation (4.3 1 and equation (4.4f I the behavior of the solutions e i 
around the corresponding singular points z± ^ ^ z± (rip ; r-|_ , r_ ) is defined by the dominant factor 
^CT0/3±/2^ All other factors in equation (4.3) are regular around these points. The same solutions are in 
general singular around the corresponding singular points z± = 1 = z±(r±;r+,r_. Hence, at this point 
we have a complete analogy with the Regge- Wheeler case, including our notations. 

Only two of the sixteen solutions (4.3 1 are linearly independent. For different purposes one can use 
different pairs of independent local solutions. 

Using the known asymptotic expansion of the confluent Heun function [ID] we obtain, as a general- 
ization of (2.51, two asymptotic solutions of Tome type. These are local solutions of TRE around its 
irregular singular point |r| = cxd in the complex plane C^: 

.±oo / \ i(TQ a; ( r+(r+ +r_ ) In r ) \ ■ / i ' 

(r; r+, r_) ~ e ^ + ' ^ > a,- | ±- 



,m.(Ta ,(7/3 ,<y-f 



j>0 



] + l + {l + <Tc)s 



a-o = 1- 



(4.5) 



The notation ±oo in (4.51 denotes the two directions: 



and 



approaching the irregular singular point |r| = oo in the complex plane 



-oo on the real r-axes for 
For the coefficients aw — 



one has a recurrence relation |10j which shows that they increase together with the 
integer j. Hence, the asymptotic series (4.5) is a divergent one. 



As seen from (4.4), 



(r; r_, r_|_). Hence, one can introduce 



a new parity property of the solutions and construct a symmetric and anti-symmetric (with respect to 
the interchange r-|_ ^ r_) solutions of TRE: 



5-K, 



SYM 
,E.m 



sR 



ASYM 



^(r;r+,r_) = 
(r;r+,r_) 



r: 



(r;r+,r_) 



si?: 



^(r;r+,r_)) 



(4.6) 



Clearly, these solutions are singular at both horizons in the general case, but when one considers the 
two-singular-point boundary problem |10j on the interval [r_, r+] in KBH interior, the solutions (4.6 1 may 
be regular at one, or at the both ends for some values of the separation constants lo and E. Since this 
boundary problem is still not studied, at present we are not able to make more definite statements about 
this case. 



5 Classification of the solutions to TRE based on the ^-condition 

For TRE the (5-condition reads: 



, ,± r ± Q 



{Mt.^^^.^^^^+is'nt,.,,.,)^ (5-1) 



where 



s. 
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The calculation of the values of the coefficients in equation (5.1 1 yields two very different cases 
1. In the ffist case from relations ± ± 



Cj — we see that one is not able to fix the frequencies 
sUj^i ± ± ± and ± ± ±- Instead, choosing aa = crp = — — sign(s) 
the degree of the polynomial Ajv+i-condition in the form 



-a and using (5.11 one fixes 



,A^+1 =2|s| 



(5.2) 



2. In the second case the coefficients iLj 



a) For £+ 



are not zero and one can fix the values of the frequencies 
from equation (5.1 1. Thus one obtains two different types of exact equidistant spectra: 



=F,±,± 



±2, M+ ±_± = M-_±_± 



(5.1) produces the pure imaginary equidistant frequencies 



and ?^^,±,± = ^=p,±,± = (-^+1) the (5-condition 



sUJ 



±i 



N+1 

AM ' 



> - integer. 



(5.3) 



Note that these frequencies do not depend on the spin-weight s and azimuthal number m, nor on the 
rotation parameter a, i.e., this spectrum is not influenced by the rotation of the waves and the rotation 



of Kerr metric. The frequencies (5.3 1 are defined only by the monopole term in multipole expansion of 
the metric. 

b) For all other cases the coefficients ^.^ , and Mj_^ are not constant and one obtains the 

following two similar double-equidistant spectra of frequencies: 



,+ 



1 I f 

^'^N,m,^,^,± = ^^w,m,T,±,T = ^ 4M i ^ ^ 7^ ' ^ ^ ^ ' ^ " ^' " mtcgcrs, (5.4a) 



=^ mn_ ± ^ -1 ] {N + l±s), N>0, m - integers. (5.4b) 



6 Polynomial Solutions to the TRE 

The (5-condition is not sufficient to ensure polynomial character of the solutions, but it yields the basic 
classification of the solutions, described in previous Section 5. In accord with it one obtains two classes 



of polynomial solutions to the TRE, imposing in addition the Ajv+i-condition (1.9) 



6.1 First Class of Polynomial Solutions to TRE: 



The solutions of this class correspond to the first case in the Section 5 and obey the equation (5.2). The 



inequality ^iV = 2|s| — 1 > excludes the existence of scalar perturbations (|s| = 0) of first polynomial 
class, just as in the case of RWE. 



6.1.1 The General Case: 

For brevity we denote the solutions ^ „ _cr -o- -cr(''' ''+' s-^w _e m(''' ^'^^ them the 
parameter /i takes the values n = slJ-^ ^ m(''+' = 1, • ■ ■ , 2|s| - the solutions of the algebraic equation 
(1.9), which now takes the form: A^^|(^) = 0. Its degree is 2|s| = 1,2,3, or 4, depending on the spin of 



the perturbations \s\ — 1/2,1,3/2,2. Making use of (1.7b I, and (4.4al-(4.4e I, we obtain for the separation 



constant E = sE^ ^ m(''+' ''-)' k=l, . . . , 2|s| the expressions 

sE^,k,mir+.r-)^st^^^k,mir+,r^) + \s\i\s\ - 1) - auj{a^ - 2m) + 2ia{2\s\ - l)ujr^, 
Applying the explicit expressions for the roots sfJ-^ k m(''+'''-)' obtain: 



sE^^^{r+,r^)^ -aLu{auj - 2m) 



for \s\ 



(6.1) 



(6.2) 
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r_) = —aLu(auj — 2m) — 2(— 1)*^ auj{au! — m) : for |s| = 1, /s = 1, 2. 



(6.3) 



For the gravitational waves (|s| — 2) one has to find the quantities sA* J ^ solving algebraic 

equation of fourth degree (/i) — 0. The explicit form of its roots is too complicated and not necessary 
for the purposes of present article. It is more instructive to demonstrate here the result, obtained using 
the Taylor series expansion of the solutions ^ „j(r+, r_) around the zero frequency lo — Q. 

Thus we obtain for \s\ =2, fc = 1, 2 the eight values: 



sE^. =2-4 



ni-i{~l)'''^\ aw+6(m^ + i(-l)''2m ( (rn^ 



1) 



M 



2M\ 
a J 



UP 7 



-j{a^f+ (6.4) 



For |s| = 2, m^O, k — 3, A we have another eight values: 

'3A/2 2\ 1 



.E^ 



i {-!)'' A^/^mlaLjl 1 + i3 1 + 



8a2 



2 

" 3 



niau! +02{aLLj) + 



/3M2 



V a2 



(6.5) 



Clearly, these series describe two kinds of solutions with a completely different behavior around the 
origin a; = 0. In particular, the series (6.4) and (6.5 1 have different limits: 2 and 0, correspondingly, when 



0. For the solutions (6.5) the origin = is a branching point, etc. 



The independence of the values of gE^ k m '^^ (|6.4| and (|6.5[) on the ± signs is a result of the polynomial 



character of the solutions, i.e. of the regularity of the corresponding HeunC-factor simultaneously on both 
horizons r±. 

For a complete solution of the problem one has to determine the frequency uj. Hence, one needs 
additional relation between the parameters E and lu. Such relation may appear when one solves the TAE 
(See next Section 8.). 



6.1.2 The Special Case of Schwarzschild metric: 

For the special value of the parameter a = we have r_ — 0, — 2M. This is the case of perturbations 
to non-rotating SBH, now described in terms of Weyl scalars. For consistence with the description of 
SBH perturbations by RWE here we have to use units in which 2M = 1 (See the Introduction.). The 
parameters in the solution (4.3 1 acquire the limiting values 



"4 



2zu;, Z?-!- = s, 7+ = s + 2iuj, (5+ = 2iLu{s — icu),?]^ = 
a_ = —2iLU, /3_ = s + 2iLU, 7_ = s, (5_ = —2iLu{s — icu), ri_ = —E - 



-E+^ 



s 



2uj + 2isuj. 



(6.6) 



These differ from the values of the parameters (2.4a I for description of the perturbations to Schwarzschild 
metric in the RWE approach. 

In the limit a — s- the equation (5.1) does not define the frequency if = = — ^cT; 
because then one obtains Lr^^j^^^^ — 0. If, in addition, a = sign(s), then the i5-condition is fulfilled 



for the special polynomial solutions of the first class, denoted as 

The equation (5.1 1 yields the relation sA^=|s| — 1>0 - precisely the same as in the case of polynomial 



.('^;i,o). 



solutions of first class to the RWE (See Section 3.0.1.). Scalar perturbations of this type do not exist. 



In case of integer spins \s\ 



( |1.7b[ ), and ( |4.4a| -( |4.4e[ ) with r+ 
E^,E^,^^, fc = l,...,|s|: 



1, 2 the roots fi 

1, r_ = and a 



u,k,m ' 



.E^ 



.E^ 



for |.s| = 1, 



= 1, . . . , |s| of the equations A|^|(^) = 0, 
produce the following simple expressions for 



(6.7) 
(6.8) 



(-l)''Vl - iGcruj : for \s\ = 2, fc = 1, 2. 
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For a complete solution of the problem, one needs an additional relation between the parameters E 
and oj. Such relation may be found solving the TAE (See Section 8.). 



6.2 Second Class of Polynomial Solutions to TRE: 



According to section 5 the solutions of this class originate from the second case of (5-condition and fall into 

■ formulae ( 5.3 1 and ( 5.4 1 , 



two subclasses: a) and b) . The complete definite frequencies lu 

yield algebraic equations A^^j(/x) = with {N + 1) roots fj. = sM^,„,„,^^,„_3,^^(?'+,r_), n = 0, 1, . 
It seems hard to derive explicit analytic expressions for these roots, but their numerical values can be 
easily obtained. Using the values of ™ 

complete definite values for the parameter E = gE^ ^ ^ o-q^ (''+>''-)■ 
a) Putting a = sign(s), in the case of the frequencies (|5.3[) we obtain 



(r_)_,r_) and equations (1.7b I, (4.4al-(4.4e I we obtain 



.E 



n,m.,a. — (T, — a 



7i.ni.a, — a. — a 



4wV^+2icrw(2Af|s| - r±). (6.9) 



,E 



b) In the case of the frequencies (5.4a I, (5.4b I we have correspondingly 
± 



N,n,m, + ,-,+ 



i 



(6.10a) 



+ 8m (1 + AfVa^) aw + (l + 10rzp/r± + 9 ir^/r±f - 4 ir^/r±f) wV^) /p^, 

± - _^-,M^„,„__^^-±*2(7n+2ac^(l-2MVa'))/p- 

2psatJ - 4:(m^ + 2m {l-3M^/a^) auo - {l-bM^/a^) {auf^ (6.10b) 



Tp± ..It 



With cj and £^ given by the equations ( |5.3[ ), ( |5.4[ ) and ( |6.9[ ), (6.101 we have no more free parameters 
in the problem at hand. As a result the corresponding solutions to TAE are fixed unambiguously by the 
designated group of equations, obtained for the second class of polynomial solutions to TRE. 



7 Exact Solutions to the Teukolsky Angular Equation in Terms 
of the Confluent Heun's functions 

In terms of the variable x ~ cos 6 the TAE has three singular points. Two of them: x_ = — 1 (i.e., 9_ ~ tt 
- "south pole") and = 1 (i.e., 0+ = - "north pole") are regular singular points. The third one 
Xoo = oo is irregular singular point. It is remarkable that introducing the notations 

z+ = z+{9) = {cos{9/2)f , z_ = z_(6') = (sin(6l/2))^ z+ + z_ = l; (7.1) 

and 

a± — ±4aa;, 6± = s =p 777, c± = s ± 777, d± ^ ±isauj, n± ~ =p 2saijj — aio^ — E. (7.2) 

we can write down the 16 local solutions of the TAE in the form 

« S^.E^ra^a,^ .a,,a^= e"" ^ ^ HeunC(aaa^ ,abb^, a^c^ ,d^,n^,zj, (7.3) 



which is very similar to the form of the (2.3 1 and (4.3). 

Note that according to the ( |7.3| the behavior of the solutions e m a a around the corresponding 



singular points z± — Q — z±{6^) is defined by the dominant factor z± ^ . All other factors in (7.3 1 are 
regular around these points. The same solutions are in general singular around the corresponding singular 
points z± = 1 = z±{6±). Hence, at this point we have a complete analogy with the cases of RWE and 
TRE, including our notations. 



15 



Only two of the sixteen solutions (7.3 1 are linearly independent. For various purposes one can use 
different pairs of independent local solutions. 

In the case of integer spin weights s — 0, ±1, ±2 there exist an additional complication. The confluent 
Heun functions HeunC(Q;, /?, 7, (5, 77, z) are not defined when /3 is a negative integer |10j . Therefore if we 
have a negative integer (3 = a{,b^ < 0, we must write down the corresponding solutions in the form 

sS^,E,m^a^,a„a^ ^ e"' ~ ~ HeunC(gaa^ , abb^ , acc^ ,d^,n^,zj. (7.4) 

For this purpose we define the concomitant confiuent Heun function 

HeunC(a,/3,7,(5,?7,z) = z"^HeunC(a,-/3,7,5,77,z) f e °'^C^ - C) '' ^ ^^ ^^ 

(HeunC(Q!, — 7, (5, 77, z)) 

Note that this function is well defined for non-positive integer /3 = abbj_ < 0, together with the confluent 
function z~''HeunC(a, — /3, 7, (5, 77, z). In this case the confluent function 2:~'^IIeunC(Q!, — /3, 7, 5, 77, z) rep- 
resents the local regular solution around the singular point z = and the concomitant confluent function 
HeunC (Q:, /3, 7, <5, 77, z) represents an independent local singular solution around this point. It can be shown 
that for negative integer /3 the concomitant confluent Heun function possess the form 

l/3| ^ 

HeunC (a,/3,7, J, 77, z) = -|- hi{z) + h2{z)ln{z), all c„ 7^ 0. (7.6) 

n=l ^ 

Here hi 2{z) denote two functions of the complex variable z, which are analytic in vicinity of the point 
z = 0. In the problem at hand \(3± \ = \s^m\. The logarithmic term presents in the concomitant confluent 



Heun function when |/3| — 0, too, but then we have no poles in the solution (7.6). Its form otherwise is 



similar to (7.6 1. One can reach the last results using general analytical methods, described, for example, 
in i29j. 

8 Regular solutions of the TAE 

Using ( |7.1[ ), (7.2 1, and (7.3 1 we write down two independent solutions to the TAE in the form: 



._/e--Hf)^(cosi)l^-"l (sini)^+'^ 



(-"5) (sin|)l^+™l(cosf)^ 



ReunC \ AaiL!, \sT m\ ,s ±m,±AsaiL!,{rn'^ + s'^)/2^2sauj - auj"^ ~ E, ^'^'^^ I. (8.1) 
Then the physically obvious symmetry of the problem becomes transparent: the solutions sS^^^^{9) 



2 ' 



and sS^ ™ (^-1 ' interchange their places after the substitution s —i- —s, 9 ^ n — 9. 

One will have solutions sSJ'%°^{9), regular at both poles, if and only if sS^^^^^{9) — const sS^ "^^^{9), 
or, equivalently, if the Wronskian vanishes: W sS^^^^j^{9), sS^^^j^{9) 
E — E(auj, 771, s) in the form 

E = (777^ 4- s^)/2 - acj^ + e{auj,m, s). 



— 0. This condition determines 
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Here e{auj, to, s) is a solution of the transcendental equation 

HeunC'(4aa;, |to + s|, s — m, —4sauj, +2sauj — e, (sin |)^) 
HeunC(4ati;, |to + s|, s — m, — 4sacj, +2sacL; — e, (sin |)^) 
HeunC'(4aw, \m — s\, s + m, +4sacj, — 2saa-' — e, (cos |)^) 
HeunC(4aCi;, |to — s|, s + to, +4saaj, — 2.saaj — £, (cos |)^) 

2(cos|)^ 2(sin|)^ 

valid simultaneously for all values of G (0, it). The symbol HeunC' denotes the derivative of the confluent 
Heun function. 

Since the Eq.(8.2) is fulfilled for tt — 9 e(0,7r), if and only if it is fulfilled for 9 €(0,7r), its symmetry: 
s—^—s is obvious. The properties of the function HeunC |10j yield the additional symmetry: ni-^—m of 
the Eq.([8^. 

Let us consider the limit acj of the regular solutions of TAE. Since 
HeunC(0, /3, 7, 0, ry, z) = (l -z)/5+7+i+\//3'+7^+i-4r, ^ 

2-f'i ^ , ^ ;p+i;2 ]\>i-'i) 



in this limit the Heun functions in (8.1 1 can be reduced to the Gauss hypergeometric ones. Thus we 
obtain the solutions in the form 



Here 



sSo,E%i^) ■= { ) . li\s+m\; g!s-^ f 2i^l ao,±:&o,±;Co,±; ) . "^{2 I • (8.4) 
Us™ 2) (cos 2) J V 



2 ~ 

(sinf)^ 



Is^ml + sim 1 r— — , |sTTO| + s±m 1 r— — , , , , 

^ + - + ^E+l/A, = ^ +--v/i?+l/4, c„,^ = |sTm| + l. (8.5) 



The condition ( |8.2[ ) reduces, correspondingly to 
2-^^! (ao,+ >^o,+ ;Co,+ ; (cosf)^) s-F'i (ao,-:^,-;Co,-; (sinf)^) |s-TO|-(s-m) |s+m|-(s+TO) 



\ 2 



= 0. 



2-Fi (ao,+ -^o.+ ;co.+ ;(cosf)') 2^1 (ao.->&o,-;Co,-;(sinf)') 2(cosf)' 2(sin|)^ 

(8.6) 

Using the well known properties of the Gauss hypergeometric function 2F1 one can derive from ( |8.6[ ) 
the spectrum E{Q\ I, s, to) = ^(Z + 1), I = l(s, to, n) = max(|TO|, \s\) + n, n = 0, 1, 2, . . . ; with real values 
of the separation constant E{0;l,s,m). The numerical analysis of the (8.2), written directly in terms of 
confluent Heun functions, conflrms this standard result for the limit auj = 0. The corresponding regular 
confluent Heun functions in ( |8.1| in this case are reduced to Jacobi polynomials. 

The solutions E{aw, s, to, I) for small aco have been studied many times |12[ 115] in the form of Taylor 
series expansion E{auj; s,m,l) = Z(/ + 1) + i?„(aa;)" without the use of the Eq.(8.2l. A little bit 
surprisingly, the solutions (8.1) with aw ^ , regular at both poles, are not polynomial and can be 
represented as an infinite series with respect to Jacobi polynomials. Further we denote these regular 
solutions by «5«f,^(0). 
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9 Classification of the solutions to TAE based on the ^-condition 



For solutions sS^ ^ (|7.3|) to TAE the ^-condition reads: 



= ^^^^+iV+l+f.„ + ^^,.. 



(9.1) 



Comparing this equation with the corresponding one for TRE - (5.1 1, we see both essential differences 



and similarities. For the coefficients in equation (9.11, which are analogous to the ones in (5.11, one 
obtains: 



<Ta ,(7b^Crc 



= 0, M±_,^_,^^,^ = T™K-a,)i =7V+l+( (7, + — ^ I s. 



1^ 



Hence: 



i) The coefficients vanish identically, in contrast to the coefficients in equation (5.11 

Consequently, there are no cases in which the condition (9.1 1 can fix the frequencies uj. 



ii) The form of the coefficients M of both equations (5.1 1 and (9.1 ) is the same only for a /M =1/2 



ill) The coefficients ?\f^ of both equations are of the same form. 
As a result we obtain two different cases: 

1. For ab = (Tc = (Ta = —u equation (9.11 fixes the degree of the polynomial AAr+i-condition in the 
same form as equation (5.2|^ 



sN+l^ 2\s\ > 1 for \s\ > 1/2. 



2. For cTf, = — (Tc we obtain 



1 = zb mcTfo — > 1. 



(9.2) 



(9.3) 



As seen, in the case of TAE the only function of the (5-condition is to relate the degree N of the 
Ajv+i-condition with the spin-weight s and the azimuthal number m. 

10 Polynomial Solutions of the TAE 

It can be shown that the polynomial solutions sS^ ^ „j ^ ^ are regular at the south pole and certainly 
singular at the north pole. For the polynomial solutions sS~ ^ ^ we have a similar result, but with 
regularity at north pole and singularity at south one. 

Studying the polynomial solutions to the TAE one must take into account one new obstacle. The 



HeunC-factors in solutions ( 7.4 1 do not become polynomials for negative integer abb^ < and for otherwise 



arbitrary values of the parameters, because of the presence of ln(z) terms in equation (7.6 1. Hence, in 



contrast to the previous cases of RWE and TRE, looking for polynomial solutions to the TAE we must 
impose the additional requirement abbj_ > in the case of integer spin weights s = 0, ±1, ±2. 



Using the relations (1.7b) and (7.2) we obtain the general formula for the constant E in the form 



=fi^ — auj^ ^ 2CTa(lT cTfc"^ + (ca + (Tb)s^aLo H — — -m {(JbTn T 1) H — — -s {abS+ 1) 



(10.1) 



Further analysis shows that we have again two classes of polynomial solutions to the TAE, as in the 
cases of RWE and TRA, but their structure in some cases may be different. 



^The alternative case crj, = = —fa leads to a non-interesting relation N + 1 = 0. 
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10.1 First Class of Polynomial Solutions to the TAE: 



,m, — cr, — (7, — o" 



with a = sign(s). For them the condition (9.2) is fulfilled 



These are the solutions 

independently of the values of the integer m, but for integer \s\ the specific requirement (7f,b^ > yields 
the restriction \m\ > \s\. As in the previous cases of first class polynomial solutions (See Sections 3 and 6.) 
the value s = is eliminated by (9.2 1. Hence, we have an infinite series of first class polynomial solutions 

Preserving the accepted in the previous sections style 
9^ = 9^ 

^ uj,E,m ^ uj^E,7n, — (T, — a, — <7' 

For them the Ajy+i-condition reads A2|s|(^) — and has 2|s|-in-number solutions 



to the TAE for all admissible values of s and m. 
we denote the polynomial solutions to TAE of the first class as 



formulae (10.1 1 one obtains 



uj,k,in' 



From the 



\s\{\s\ - 1) - auj^auj - 2m) =F 2(t(2|s| - l)aw, 



(10.2) 



where k — 1, . . . , 2\s 
Solving the An+ 



, s — ±1/2, ±1, ±3/2, ±2 and for integer \s\ in addition |m| > |s| 
i-condition, we obtain for the different values of |s| as follows: 



± 



-au^aui — 2m) 



1 
4 



for \s\ 



(10.3) 



,k,rn 



-aLo{auj — 2m) — 2{—l)^^/aLu{au; — m) : for |s| = 1, \m\ > 1, k — 1,2. 



(10.4) 



Note that the values ( 10.3 1 and ( 10.4) of the separation constant E obtained for the first class polyno- 



mial solutions to TAE are the same as the corresponding values ( 6.2 ) and ( 6.3 ) for the first class polynomial 



solutions to TRE. The author's attention to this fact was drawn by Dr. Roumen Borissov during the re- 
cent discussions on the applications of Heun's functions to TME. The important consequences will be 
considered in a separate article [30] . 

For the gravitational waves (|s| = 2) the quantities fc m solutions of the algebraic equations of 
fourth degree AJ(/x) = 0. We do not need here the exact form of these roots. It is quite complicated. 
Below we present only the form of the separation constant E for TAE, obtained making use of the Taylor 
series expansions of the roots around the point au! = 0. 

Thus we obtain for |s| =2 and fc = 1,2 the following eight values: 



.E^ 



,k,m 



2 - Amauj - i(-l)''12^(m - l)m{m + 1) (aw)^/^ + 6 ^m^ - ^ {auj)^ + 05/2{aio), 



(10.5) 



and for |s| = 2, to 7^ 0, fc = 3, 4 another eight values: 



3m 

V TO 



auj + 02{auj)^+8mauj-6(m^-^{au;)^ + 03{auj). (10.6) 



As seen, for gravitational waves of first polynomial class the values (10.5) and (10.6 1 of the corre- 
sponding constants E differ substantially from the analogous values (6.4) and (6.5) of the constants E, 
obtained for TRE in Section 6.1.1. This is in sharp contrast to the case of neutrino waves (|s| — 1/2) of 
first polynomial class and to the case of electromagnetic waves (|s| = 1) of this kind. 

It can be shown that this phenomenon reflects the difference between the Starobinsky constants for 
solutions with spin \s\ = 2 to TAE and TRE [TH|S]. For spins \s\ = 1/2 and 1 the Starobinsky constants 
for solutions to TAE and TRE are the same. 

Despite the pointed difference, the first-polynomial-class-solutions to TAE and TRE with spin |s| = 2 
have similar qualitative properties, discussed at the end of Section 6.1.1. 
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10.2 Second Class of Polynomial Solutions to TAE: 



We have a finite number of second class polynomial solutions to the TAE. For them the conditions > 
and (Tbb^ > must be satisfied simultaneously. Altogether there exist 24 such solutions gS^ e m =f ± =f ■ 



s = +2, m = -l, 0, 1, 2; s = +l, to= 0, 1; 

s = -2, m=:-2,-l,0, 1; s = -l, m=-l,0; 

s = +2, m==-2,-l, 0, 1; s = +l, m=-l, 0; 

s = -2, TO = -1, 0, 1, 2; s==-l, m= 0, 1. 



(10.7) 



The relation between the constants E and lu follows from (10.1 1, when /i in it is replaced by the solutions 
of the Ajv_|_i-condition in the form A^^_^^^^(fi) = 0. 



11 The 256 Classes of Exact Solutions to the Teukolsky Master 
Equation 

Combining the studied in the previous Sections solutions to the TRE and TAE wc can construct the 
following 256 classes of exact solutions to the TME 



s-^i^,E,m.i7a 



(r;r+,r„)s5; 



± 

uj,E,7n,cra 



..M (11-1) 



For specific physical problems one has to impose specific additional conditions, like stability conditions, 
boundary conditions, casuality conditions, specific fixing of the in-out properties, regularity conditions 
etc. Thus one selects some specific combinations of solutions to TRE and TAE in ( |11.1[ ) and derives the 
spectrum of the separation constants lj and E in the given problem. 

For example, choosing solutions to TRE which enter both event horizon and 3D-space infinity we are 
studying Kerr black holes. If in addition we choose regular solution to TAE, we will obtain QNM of KBH. 
The choice of one-way polynomial solutions to the TAE in combination with BH boundary conditions for 
TRE will produce jets from KBH |151|25]. The combination of outgoing one-way polynomial solutions of 
TRE and regular solutions of TAE seem to be proper for description of supernovae explosions and may 
produce the structure of their outbursts. The use of one-way polynomial solutions both for TRE and 
TAE [30] seems to be most natural for description of jets from object, different from KBH. Constructing 
simple models of different kinds of compact rotating relativistic objects one can use Dirichlet's, or semi- 
Dirichlet's boundary conditions on proper surface outside the event horizon, or even outside the ergosphere 
of Kerr metric, i.e. acting by analogy with the construction of such models in Schwarzschild space-time 
[3]. Excluding the ergoregion of Kerr metric from physical consideration we may ensure the stability of 
the corresponding objects, which otherwise may be problematic in some domain of parameters |14l . The 
exact polynomial solutions to the TRE of equidistant spectrum may be useful for quantum gravity, etc. 

The solutions ( 11. 1[ ) do not necessarily have a direct physical meaning. Instead, proper linear combi- 
nation of specific solutions, which obey the corresponding boundary conditions, is to describe the Nature. 
In general the solutions ( |11.1| have to be considered as an auxiliary mathematical objects - (maybe sin- 
gular) kernels of integral representations (1.3 1 of the physical solutions. The choice of the corresponding 
amplitudes sAuj.E,m,aa,<yft-iy~,-iya,,ai,,i7c will fix completely the physical model and can ensure the convergence 
of the integrals and discrete sums to physically acceptable solutions. 



12 Conclusion 

We have demonstrated that the confluent Heun's functions are the adequate and natural tool for unified 
description of the linear perturbations to the gravitational field of Schwarzschild and Kerr metrics outside 
the corresponding horizons, as well as in the interior domains. These functions give us an effective tool for 
exact mathematical treatment of different boundary problems and corresponding physical phenomena. 
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Large classes of exact solutions of all possible types to the perturbation equations, both of Schwarzschild 
and Kerr metrics, are described and classified uniformly in terms of confluent Heun functions and confluent 
Heun polynomials. Using these functions we have re-derived the known polynomial solutions and found 
a large number of new ones. 

We have to stress especially the newly obtained singular polynomial solutions of Teukolsky angular 
equation. These solutions may describe in the most natural way the collimation of the observed relativistic 
jets, related with different kind of astrophysical objects. 

The solutions of the remaining basic mathematical problems and specific physical applications of the 
obtained results will be published elsewhere. 
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